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OPTIMUM DESIGN OF MAGNETIC BRAKING COILS WITH SPECIAL 
APPLICATION TO LEWIS DROP TOWER EXPERIMENTS 

by H. G. Kosmahl 
Lewis Research Center 

SUMMARY 

General analytical treatment of a design of magnetic braking coils and magnetic brak- 
ing action is given. However, numerical results are computed and tabulated only for the 
specific application and requirements of the Lewis Research Center (400-ft) drop facility. 

A method for calculating the mutual inductance of unconventionally wound coils has 
been developed. Braking action has been computed for the case of a primary coil used to 
stop the motion of a secondary coil moving along the common axis and with a given mass 
and initial kinetic energy. The energy required in the primary coil can be minimized. 

The amount of deceleration can be adjusted by varying only the primary voltage (current). 
The preceding tasks may be accomplished with a primary system requiring a minimum 
number of turns. 

For the specific case of a freely falling payload in the Lewis 400-foot drop facility, 
weighing 1000 kilograms and having a final velocity of 50 meters per second, a primary 
power of between 10 to 20 megawatts associated with a primary current of about 100 kilo- 
amperes is required for 10 seconds to stop the payload with a deceleration equal to or less 
than 30 g's. The corresponding primary magnetic field is of the order of 1 Weber per 
square meter (10 000 gauss). 

A primary coil with a linearly increasing winding density has been found to require 
the least energy for a given maximum deceleration. Because of the nonlinear nature of 
the equations involved, additional computations may be necessary for numerical cases 
not tabulated in this paper. For the same reason, a closed form optimization procedure 
is not feasible. 


INTRODUCTION 

The use of magnetic brakes in many industrial and scientific applications is a well 
established practice. The mechanism of operation can be best explained from Lenz’s 
law which says that magnetic fields induced in one circuit as a result of interaction with 
another circuit oppose each other. 

The purpose of the present investigation is to find a configuration of a primary mag- 
netic coil (in the following coil 1) and of a secondary coil 2 such that the motion of coil 2 
through fields of coil 1 can be brought to a stop under controlled conditions requiring 


minimum primary energy. For the particular application being considered herein (i. e. , 
the drop experiments in Lewis Research Center zero-gravity drop tower facility) the 
payload falling freely through a height of 120 meters approaches a speed of about 50 meters 
per second. It must be decelerated to almost a perfect stop with the deceleration not to 
exceed a specified maximum. The solution of this and other similar problems requires 
the knowledge of the mutual inductance M between two unconventionally wound circuits 
as functions of distance between them. (Symbols are defined in appendix A. ) The deter- 
mination of M (appendix B) was accomplished by expanding a number of winding density 
functions (winding schemes) in a complex Fourier integral and applying them to a Bessel 
integral expansion of the vector potential. From the latter the mutual inductance and the 
self inductance can be determined. 

The circuit equations as well as the equation of motion of the payload in magnetic and 
gravitational fields are then solved on a computer and an optimization of the coil system 
with respect to minimum energy and minimum size is carried out. 

The significance of these results lies in the possibility of designing a system requir- 
ing a minimum of primary energy, in which the braking action is accomplished in a pre- 
cisely adjustable manner with respect to the desired degree of deceleration by simply 
varying the voltage of the primary circuit and/or the length of the primary coil. 

ANALYSIS OF MAGNETIC DECELERATION 

The interaction of two, predominantally magnetically coupled coils can be best 
described by their mutual inductance M. Figure 1 shows two magnetic coils denoted by 
subscripts 1 and 2. Coil 1 is connected to a power source with the load voltage V o and 
an impedance Rg « R^. The origin of the coordinate system z = 0 represents the 
(physical) end of coil 1, the latter being immobile. Coil 2 is short circuited and is 
assumed to move only in the z-direction, which points in the direction of the earth*s 
gravitation. 

A magnetically noninteracting mass m Q is affixed rigidly to the turns of coil 2. The 
latter has a mass mg so that m = m Q + mg represents the full mass of the second coil. 
Suppose that a steady-state current i ^ q flows in coil 1 prior to any interaction with 
coil 2. If the latter is permitted to fall in the negative z-direction due to gravitational 
force and its momentary distance from z = 0 is z, the following equations are applicable 
to determine the instantaneous velocity z|z(t)J of the second coil: 

V =l.(.)E. t L 1 ^ + «i 2t M(^ (1) 

° 1 1 1 dt dt 2 dt 
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Figure 1. - Schematic arrangement of coils 1 and 2 for falling payload. 


O = i 2 (t)R 2 + L 2 — + i (t) + M(z) ^ 

^ * * dt dt 1 dt 


( 2 ) 


The repelling magnetic force F m acting upon the turns of the falling payload is 


F m = ^ 


(3) 


where z is the distance between defined points in coils 1 and 2. The equation of motion 
of the payload is 


z 


1 1 ( t ) 1 2 ( t ) 9M(z) 
m dz 


(4) 
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Normalized force, G; secondary coil current, i 2 , kA; or velocity, z, m/sec 


Multiplication by z and integration between t = 0 and t = t (fig. 1) result in the energy 
equation 


/ t=t 

i ± (t) - i 2 (t) 

) 


3M(z) 


z(t)dt 


since z(0), the velocity at z = h, is zero. The initial value conditions are i^(0) = i^ 
and i o (0) = 0. After M = M(z) has been computed with the help of relations derived in 



Distance, x, m 


Figure 2. - Typical plots of velocity, secondary current, and normalized 
force against distance between coils. Winding density, wj = 8^1 


appendix B, equations (1), (2), and (5) 

are solved simultaneously for z(t) and 

z(t). The values z(t ) = 0 and z(t ) 

s s 

give the stopping time t and the 

s 

stopping position z(t ) = x„ of the 

s s 

payload. The total output of computer 
data contained the currents i^ and 
ig, M[z(t)] and 3M[z(t)]/3z(t), the 
normalized force G = (i^i 2 3M/3z)/mg, 
the inductances and L 2 , the re- 
sistances Rj and R 2 that have been 
computed as functions of the coil diam- 
eters 2a and 2b, wire diameters d^ 
and d 2 , the winding density w^ and 
w 2 , the coil lengths and Z 2 , and 
the specific resistivities and p 2 
of the coil wires. 

Notice that the decelerating force 
persists even after z has reached 
zero (fig. 2) because the time constant 
t 2 = L 2 /R 2 prevents an instantaneous 
decay of i 2 to zero. Thus, F m 
remains finite over a period of a few 
constants, and a ^perfect" stop is 
possible. 

A payload m Q of 900 kilograms 
has been assumed to which the mass of 
aluminum turns affixed to the payload 
has been added. Because of specific 
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requirements the diameter of the primary coil has been kept constant at 2a = 3 meters. 

The primary turns consisted of 0. 1- meter-diameter copper wire throughout with a spec- 
ific resistivity of 1. 8xlO -8 (f2)(m). 

DISCUSSION OF RESULTS 

The more important computer results have been presented in table I. In addition to 
geometrical parameters they contain the primary coil voltage Vj, the steady-state pri- 
mary current q at the beginning of interaction, the power P, the normalized force 
G in multiples of g, and the position x at which the center of the payload turns ap- 
proaches zero velocity, measured from the bottom end of the primary coil z = 0. Also 
listed are two values of the specific resistivity Pg for aluminum which have been assumed 
to be 3xl0 -8 (£2)(m) at room temperature T = 300° K and 0. 6x10 8 (S2)(m) at T = 80° K 
(liquid nitrogen (unpublished data, obtained by John C. Fakan of Lewis, show values as 
low as 0. 3xi0~ 8 (S7)(m)). It should be noted that the primary current ij is not absolutely 
constant; however, because of the relatively large time constant r^, the current i^ 
changes by less than 10 percent of the initial i^ q value. 

The maximum temperature rise in the secondary turns was less than 25° C, with an 
average of approximately 10° C, rendering this effect entirely unimportant. 

Table I refers to three different winding schemes. For deciding which scheme is 
preferable over the others the primary power must not be the sole criterion because of 
the presence of restraining conditions and ease of construction. Two important conditions 
in drop tower applications are: the time constant to be as small as possible in order 
for the current i^ to build up as close to its final value i^ q as possible in a time 
t < 10 seconds. The latter number is approximately twice the free fall time through a 
height h equal to 120 meters, the time t = 10 seconds being the maximum period avail- 
able for the buildup of current in case of a two way "up and down" experiment. 

The other conditions for the Lewis Research Center drop tower experiments is the 
upper limit of deceleration G to be smaller than 30 g*s. The computations show clearly 
that the power requirements to bring about a full stop could be reduced substantially if 
the condition G < 30 were disregarded because stopping action could then be accom- 
plished with a shorter primary coil. 

An evaluation of the tables leads to the following important conclusions: 

(a) Application of liquid -nitrogen cooling to the secondary turns reduces sub- 
stantially (50 to 70 percent) primary power requirements in the case of d 2 = 0. 05 meter 
and only slightly (about 15 percent) in the case of dg = 0. 1 meter for the assumed value of 
specific resistivity. 

(b) Primary power requirement decreases rapidly as b is increased (while 
keeping a = 1. 5 m constant) from 0. 75 to 1.0 meter. 


5 



TABLE I. - COMPUTER RESULTS AND COIL PARAMETERS 


-8 —i 

Primary coil radius, 1.5m; primary coil specific resistivity, 1.8x10" (f2)(m); primary coil conductor diameter, 0.1m. 
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(c) For equal values p 2 , d 2 , Z 2 , w 2 , and b, respectively, while keeping G < 30, 
the power requirements of the three tabulated winding schemes necessary to bring about 
a full stop are about the same. (The energy requirements, however, are not equal. ) 

(d) The winding scheme Wj = const [1 - (z/Z)] produces the shortest time constant 
= 2. 3 seconds and requires the smallest number of turns to bring about a "full stop" 

of the payload. 

The meaning of full stop, as used herein, refers to a motion in which the velocity 
decays directly to zero, without oscillations in the sign of the deceleration (i. e. , no ac- 
celeration). Notice that the position of stop x in all symmetrically wound coils lies 

s 

above the coil center z = 1^/2. This happens because in symmetric coils the mutual 
inductance M has a maximum and 3M/9z = 0 at the center of the coil. Thus 9M/9z 
reverses its sign in going through the center. Now, in the expression for the force 
(eq. (3)), ij(t) remains almost constant and the sign of i 2 depends mainly on that of 
9M/9z. Because the time constant t 2 is rarely below 0. 1 second in this particular de- 
sign, a reversal in the sign of 9M/9z in going through the center would precede that 
occurring in ig. Therefore, if the payload moved over the center, a short period of ac- 
celeration would follow: the payload velocity would increase temporarily before final 
deceleration could occur to cause the final stop. Such a "two phase" deceleration with a 
short acceleration period in between has been examined on the computer and seems to 
produce a workable mode of operation requiring approximately 10 percent less power 
than "one phase" stopping. The evaluation of the "two phase" cases, however, is more 
ambiguous than that of "one phase" stopping and has not been considered in this discus- 
sion. 


The effect of changing 1 2 on the power requirements indicates an optimum length 
1 2 < ly This result is easily explained by the fact that two competing factors produce 
a minimum power requirement: increasing l 2 increases the flux through the secondary 
and, therefore, the braking force. On the other hand the mass m 2 increases propor- 
tionally with l g and reduces the effect of the decelerating force in front of the integral 
(eq. (5)). 

The highest computed value l 2 is J 2 = 3. 6 meters at which number the primary 
power requirement is still decreasing slightly; however, the dependence of P against 
1 2 flattens out and an increase of power requirements with increasing Z 2 is likely to 
occur at some large value of Z 2 , which is of little practical interest to the desired appli- 
cation. 

A clear dip in power requirements appears if w 2 is changed, as can be seen from 
table I. This optimum occurs because of the competing effects between the secondary 
time constant r 2 which prevents a rapid build-up of i 2 and of the mutual inductance 
Mj 2 . Since both t 2 as well as Mj 2 increase with increasing w 2 , a minimum in 
the required power occurs if the total interaction time is not substantially larger than r 2 . 
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Figure!. - Normalized force against distance for three different winding schemes. Secondary coil wind- 
ing density, 6 turns per meter; secondary coil length, 0.6 meter; secondary coil radius, 0.75 meter. 


Because of the approximate constancy of i 1 the primary energy requirement is pro- 
2 

portional to i^T^, where is the braking time. The smallest braking time is accom- 
plished with a steplike- shaped deceleration of maximum permissible amount over the time 
Tjj and which is zero elsewhere. Figure 3 shows the normalized magnetic deceleration 
G against distance for two different geometries. The length of coil 1 wound according to 
w^ = 10(1 - z/jj) is l j = 8 meter; that of coil 2 wound according to 
Wj = 2(1 + 5 cos ffz/ij) is ij = 12 meter. The difference in coil length is necessary in 
order not to exceed G = 30, as may be seen from comparing curves 2 and 3 in figure 3. 

A comparison of the curves shows that curve 1 except for the exponentiallike tail, more 
closely approaches the steplike function than curves 2 and 3. The latter resemble rather 
resonancelike curves. 




Figure 4. - Derivative of mutual inductance against distance for four different winding 
schemes. Secondary coil winding density, 6 turns per meter; secondary coil length, 
0. 6 meter; secondary coil radius, 0. 75 meter. 


Figure 4 shows 3M/3x against x for a variety of winding schemes; x is the dis- 
tance from the center of coil 2 to the center of coil 1 in symmetrically wound primary 
coils and from center of coil 2 to the remote end of coil 1 for the winding scheme 
w^ = 10 [1 - (z/Z j)]. A uniformly wound solenoid produces the narrowest, the 1 - (z/ Zj) 
winding scheme the closest to rectangular or steplike shape. Figure 4 shows that 
3M/0x is much wider than the deceleration G in figure 3. Both are plotted on the same 
scale of x. This behavior results from the existence of a time constant Tg which is not 
substantially smaller than the braking time t^. Thus, there is a delay before the maxi- 
mum of ig can develop and the width of G is shortened (fig. 2). 

Figure 5 shows Pp Tp G may and the primary number of ampere turns plotted 
as functions of w^ q for the winding density w^ = w^ q[1 - (z/Z^)]. As w^ q increases 
the power necessary to produce a full stop, Pj decreases over the examined range. On 
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45 


Coil 1 time 
constant 



Steady-state winding density, W]_ g, turns/m 

Figured - Power, time constant, normalized force, and number of ampere turns 
as functions of steady-state winding density for primary winding density required 
to stop payload. Primary coil length, 8 meters. Secondary coil parameters: 
winding density, 6 turns per meter; coil length 1.2 meters; coil conductor 
diameter, 0. 1 meter; specific resistivity, 3x10“° (fl)(m). 


the other hand 
turns 


G remains constant (G = 30g) as does the total number of ampere 


i 



It can, therefore, be concluded that as long as ^ and the total number of ampere 

turns are kept constant G mQ ,„. will not change. 

max 

A practical limitation in the reduction of P 1 is twofold: the increase in r 1 (which 
could be irrelevant in other cases) and the necessity to increase the radius of the turns. 
Also, at some point the change in i^ must be considered, and the integral 


f VjijftJdt 
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must be used for primary energy calculations. The power curve in figure 5 flattens 
out. Its minimum lies beyond the range which is of interest for the considered ap- 
plications. 

Figure 2 (p. 4) is a plot of the velocity z, the secondary current ig, and the force 

G against x. Note that at x = x c the payload comes to a stop (i. e. , z = 0), but both 

s 

ig and G are different from zero, which makes a perfect stop possible. The shape of 
G follows closely that of ig. Because of the constancy of i^ this is only possible if 
9M/0x remains approximately constant, in agreement with curve 4 in figure 4. 


CONCLUDING REMARKS 

A design of a magnetic coil system in applications as a magnetic brake has been in- 
vestigated and a method for calculating the mutual inductance and self -inductance of non- 
conventionally wound coils has been developed. The main features of this brake are 
minimum primary energy requirement to bring to a stop a payload with a given initial 
kinetic energy, production of a smooth deceleration in a steplike function with the degree 
of deceleration being easily adjustable by changing only the primary voltage and/or pri - 
mary coil length Ip accomplishment of the preceding tasks with a primary coil having 
as few turns as possible, and a required flux density equal to about 1 Weber per meter 
squared (10 000 gauss). 

A very definite optimum design has been determined under the imposed restrictive 
conditions (i. e. , maximum deceleration and coil diameter ratio). The choice of param- 
eters for minimum power requirement is generally noncritical; however, deviations of 
a factor of 2 in some of the parameters from their optimum value may render the system 
infeasible. 

It is concluded that an elegant, efficient, and reliable system could be developed 
unless stability studies, not carried out in the frame of this investigation, produced in- 
formation to the contrary. 

Although not discussed herein, it is conceivable that an electromagnetic system could 
be designed to fulfill the double purpose of launching and stopping the payload within the 
very same coil, thus eliminating the need for mechanical launching guns. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 9, 1965. 
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APPENDIX A 


SYMBOLS 


A <p 

azimuthal component of vector 

l 

potential £ 

M 

a 

radius of primary coil, m 


B 

2 

magnetic flux density, Wb/m 

m 

b 

radius of secondary coil, m 

N 

C, D 

general constants 

P 

C(k) 

wave number dependent pa- 
rameter 

R 

r 

d 

distance between coils, m 

S(k) 

dS 

surface element 

t 

ds 

line element 

V 

F 

m 

magnetic force 

w 

G 

normalized magnetic force, 
F/mg 

X 

g 

H 

9. 81 m/sec 2 

magnetic field strength, A/m 

a Os), $(k) 
6 
P 

P 

h 

height, m 

modified Bessel function of 

■fVi 

first kind n order and 


argument kr 

^0 

i 

current, A 

T 

j 

current per unit length, A/m 

T b 

K^kr) 

modified Bessel function of 
second kind, n 11 order and 

Subscripts: 

e 


argument kr 

i 

max 

k 

L 

Fourier integral wave number 
self inductance, H 


length of coil, m 

mutual inductance between 
coils 1 and 2, H 

mass, kg 
number of turns 
power, W 
resistance, 
radius, m 

Fourier amplitude function 
time, sec 
voltage, V 

winding density, turns /m 

distance, m 

axial coordinate 

phase angle 

Dirac’s delta function 

specific resistivity, (0)(m) 

permeability 

47TX10" 7 H/m 

time constant, L/R 

breaking time 

external 

internal 

maximum 
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s stopping 
cp azimuthal 

0 initial 

1 primary coil (coil 1) 

2 secondary coil (coil 2) 


Superscripts: 

vector 

first derivative 
" second derivative 
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APPENDIX B 


DERIVATION OF EXPRESSIONS FOR MUTUAL INDUCTANCE AND 
SELF-INDUCTANCE OF COILS WITH CYLINDRICAL SYMMETRY 

If it is assumed that the current sheet has a negligible thickness and circulates only 
in the azimuthal direction <p at r = r Q , only one component of the vector potential 
will exist. This can be accomplished by applying a suitable winding technique or current 
feeding technique. 

For a coil of radius r Q the latter may be represented by references 1 and 2. 

i[kz+o;(k)] 

C i (k)I 1 (kr)e dk for r < r Q (Bl) 


i[kz+Q!(k)J 

C e (k)K 1 (kr)e dk for r > r Q (B2) 

The boundary conditions for B and H on the interface r = r Q between the internal 
and external region, denoted by subscripts i and e, respectively, are 



<P: 


i = \,,e< B i = B e> 


(B3) 


and 


— — — (rA A-— — — (rA \ = j (z) at r=r 

n. r 9 r rV) Mo rarV <?> e J W' 


(B4) 


where j (z) is the surface current density per unit length 
s, (p 


h, (p^ = f ( z ) 6 ( r “ r Q ) for 0 < z < l 


and j (z) = 0 elsewhere. If only the winding density w(z) but not the current > in the 
s j cp 

turns is a function of z, the following expansion of j is appropriate: 
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V = > w < z > 


= 



S(k)e ikz dk =ywj 



|S(k) |e i [( kz+$ ( k ).l dk 


(B5) 


S(k)=± e- ikz w(z)dz= |S(k)|e i<& < k > 

27 T J 

OO 


(B6) 


expressions (B4) to (B6) the following equation is obtained for p. = M e = 


jc i(k) ± [rl j(kr)] 


- C e (k) ± [rKjOtr)] 
r 


e i[kz +a (k) ] dk 


V. 


-/ 


r ~ r o 


- CO 

Li 


C i (k)kI Q (kr o ) 


c e (k) 


kK 0 (kr 


°)] 


e i[kz+a(k)] 


dk 


= ^o V (Z) 


JUW 



|S(k) |e i t kz+ ®( k )] dk 


= 2 ^o w l J 


/ OO 

|S(k) | cos jjtZ + $(k)J 


dk 


(B7) 
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because the integration over sin from -°° to °° is zero. Equations (B3) and (B4) can be 
satisfied only if 


Qf(k) = <t>(k) 

(B8) 

- C e (k)K 1 (kr o ) = 0 

(B9) 

C i (k)I 0 (kr Q ) + C e (k)K 0 (kr 0 ) - |S(k) | 

k 

(BIO) 

With 


^(krjKjtkr) + I 1 (kr)K 0 (kr) ^ -1 

kr 


the constants and C g are determined from (B9) and (BIO) to 


c i( k ) = ^o/WilS^koKi^o) 

(Bll) 

c e ( k ) = Mq/Wj |S(k) |r Q I 1 (kr o ) 

(B12) 


Thus, the expressions for the vector potential (Bl) and (B2) may be written 

|S(k) |Kj(kr Q )Ij(kr)cos[kz + $(k)] dk (B13) 


|S(k) |lj(kr 0 )Kj(kr)cos[kz + $(k)] dk (B14) 


Mutual Inductance Between Two Coils 

Now the relation (B13) and (B14) is applied to calculate the mutual inductance M 
between coils 1 and 2, each having a constant diameter 2a and 2b, respectively. In 
addition, the winding density on the second coil is assumed constant (e. g. , Wg = const) 


V, i = 2 »o r o Jw l [ 

•'fl 


and, similarly, 


e = 2 “ 


o r c/ w l / 

Jr\ 
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even though this restriction is not, generally, necessary. 

In a distance between z and z + dz, counted from a fixed point z = 0 in the primary 
coil (subscript 1), there will be w 2 dz turns of secondary coil (subscript 2). Thus, the 
mutual inductance between the element w 2 dz and the entire flux coming out of the pri- 
mary coil is 


dM = w 2 dz 



5 l( z W>* “ dS 2 = w 2 dz 



X l< z >r=b 



(B15) 


where S 2 and s 2 designate the area and the length element of turns in the secondary 
coil, respectively. Because A^(z) is constant over the azimuth at a given radius, the 
following may be written for dM 


dM 


27Tb A i, ^,( z ) r= b w 2 dz 


(B 16 ) 


The integration (B16) over the entire length 1 2 results in 

•d+(Z 2 /2) 


M = 4/u. o abw^w 2 


p oo /»CH 

7T / dk / * |S(k)|K 1 (ka)I 1 (kb)cos[kz + $(kjjdz (B17) 

Jk=0 «/z=d-(l 2 /2) 


where d is the distance between a fixed point in coil 1 and the center of coil 2. 


Self-Inductance of Coil 

The relations developed for calculating M will now be applied to determine the self- 
inductances of axially symmetrical coils of constant diameter. This is done by first find- 
ing the flux linking any turn and then integrating over all turns. 

Modifying equation (B15) to apply to the very same coil yields 
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L = / w l< z > dz [ B l< z )r=a ,E l 

"'ll JS 1 

m / W!<^y Aj(z) r=a 

= 2vra / w 1 (z)A 1 


(z) dz 
<p v 'r=a 


(B18) 


Fourier Integral Expansion of w^z) 


The function 


w x (z) = f(z)6(r - r Q ) 


for either 


and 


f(z) ^ 0 in 0 < z < U , or f(z) ^ 0 in - 1 < z <- 

L 2 2 


f(z) = 0 elsewhere 

will be expanded in a Fourier integral: 

f(z) = f S(k)e i]sz dk 

- OO 

Note in the following that the center of expansion z = 0 coincides with the center of coil 1 
for symmetrically wound coils and with the bottom of coil 1 in all other cases, in agree- 
ment with figure 1. The distance x is always counted from the proper center of expan- 
sion on coil 1 to the center of coil 2. 
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lllilll! I II II! 


and 


w(z) 


=«i [ 1 + D (i-f M)] 

= 0 elsewhere 
k l 


- 1 < z < ^ 
2 2 


D(k) = 


sin — - 

1 A 

is k \ 


cot^l 


$(k) = 0 


= w , (l - - \ for 0 < z < 1 

°> 1 V i) 


w(z) 

= 0 elsewhere 


S(k) = |s(k)|e i$(k) 


|q/b-\ I L_ I~ 1 - cos k l _ sin k l n 

I()| j2irk[ (kl) 2 kl 2 J 


. 1/2 


tan <&(k) = 


sin kl - k l 
1 - cos kl 


Formulas for M and L for the analyzed function w 1 (z) are as follows: 

w^(z) = const 
W£ = const throughout 


Mi 


l 


klj k I 2 


dk 


12 = 8ju. o abWiW 2 I Ii(kb)Ki(ka) sin sin cos kx — 

- I 22 ,2 


where 


a > b; x is distance between centers of 1 and 2 


(B21) 


(B 22 ) 


(B23) 
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(B24) 


where 


a > b; x is distance between centers of 1 and 2 




k Z « kl q die 

x I, (kb)K., (ka)sin — - sin — — cos kx — 
1 1 2 2 k 2 


where 


(B25) 


(B26) 


a > b 
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= 8M Q a 2 Wj 


f [ 1+ — 

/ k ^1 
(CSC 

k Z A 
- cot — ± 

1 L w i 

\ 2 

2 / 


I 1 (ka)Kj(ka)sin 2 — - — 


w < z ) = w i,o( 1 -f-) 


for 0 < z < Z 


1 


M 


J f°° k l 

I f(kZj)I 1 (kb)K 1 (ka)sin — - cosj^kx + $(k)J — 

0 k 


where 


a > b; x is distance between x = 0 on coil 1 and center of coil 2 


L 1 = 8jU 


2 2 f 00 f(kZ.) kZj 

o a w l, 0 / — ~ I 1 ( ka )K 1 (ka)sin — 

Jo k 


cos 


'kZ« 

— - + $(k) 
2 


dk 


i fl a 2 Wj o / — — — I 1 (ka)K 1 (ka)sinj^kZ 1 + $(k)Jdk 

Jo k 

o O f°° f(kL) kZ, |"kl 1 

_a z w7 n / — I^kajK^kaJsin — - sin — ± + 

’ J k Z : 1 1 2 2 


dk 

,.2 


i 1 - cos k Z- sin kl 1 
f(k L) = — — — + 1 1 

1 / — O 


VU 2 (M«) 


2 kl. 


1/2 


sin kZn - kl 1 

tan <&(k Z-.) = for - — < <&< 0 

1 1 - cos kZj 2 


(B27) 


(B28) 


(B29) 


(B30) 


(B31) 
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